Rician Noise Removal by Non-Local Means Filtering for Low Signal-to-Noise Ratio MRI: Applications to DT-MRI

Diffusion-Weighted MRI (DW-MRI) is subject to random noise yielding measures that are different from their real values, and thus biasing the subsequently estimated tensors. The Non-Local Means (NLMeans) filter has recently been proposed to denoise MRI with high signal-to-noise ratio (SNR). This filter has been shown to allow the best restoration of image intensities for the estimation of diffusion tensors (DT) compared to state-of-the-art methods. However, for DW-MR images with high b-values (and thus low SNR), the noise, which is strictly Rician-distributed, can no longer be approximated as additive white Gaussian, as implicitly assumed in the classical formulation of the NLMeans. High b-values are typically used in high angular resolution diffusion imaging (HARDI) or q-space imaging (QSI), for which an optimal restoration is critical. In this paper, we propose to adapt the NLMeans filter to Rician noise corrupted data. Validation is performed on synthetic data and on real data for both conventional MR images and DT images. Our adaptation outperforms the original NLMeans filter in terms of peak-signal-to-noise ratio (PSNR) for DW-MRI.